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In order to understand galaxy evolution through cosmic times it is critical to derive the properties of the molecular gas con- 
tent of galaxies , i.e. the material out of which stars ultimately form. The last decade has seen rapid progress in this area, 
with the detection of massive molecular gas reservoirs at high redshifts in submillimeter-selected galaxies and quasars. In 
the latter case, molecular gas reservoirs have been quantified out to redshifts z>6, i.e. towards the end of cosmic reion- 
ization when the universe was less than one Gyr old. The recent discovery of molecular gas in more normal galaxies have 
extended these studies from the most extreme objects in the unviverse (SFR~1000Mq yr _1 ; quasars and submillimeter 
galaxies) to more 'normal' starforming systems at redshifts 1.5-2.5 (with SFR~1OOM0 yr _1 ). However, detecting the 
molecular gas reservoirs of high-redshift galaxies that only have moderate star formation rates (~10Mq yr~\ similar to 
the faint galaxies seen in the Hubble Ultra Deep Field) will likely have to await the completion of ALMA. 
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1 Recent Progress in understanding Galaxy 
Formation 

The last decade has seen dramatic advances in our under- 
standing of galaxy formation. Cosmic geometry, the mass- 
energy content of the Universe, and the initial density fluc- 
tuation spectrum, have been constrained to high accuracy 
(e.g., Spergel et al. 2007). Given the underlying distribu- 
tion of dark matter, the main challenge of galaxy forma- 
tion studies is to understand how the baryons assemble into 
the observed Universe. In this context, structure formation 
through gravitational instabilities has been calculated in de- 
tail through numerical studies, and observationally verified 
through studies of galaxy distributions (e.g., Springel et al. 
2006). The history of star formation, and the build up of 
stellar mass, as a function of galaxy type and mass, are now 
constrained out to the tail-end of cosmic reionization, within 
1 Gyr of the Big Bang, and pushing toward first light in the 
Universe (e.g., Bouwens et al. 2010). 

The principle results of these studies can be summarized 
as follows. The comoving cosmic star formation rate density 
increases by more than an order of magnitude from z ~ 
to 1, peaks around z ~ 2 to 3, and likely drops gradually 
out to z ~ 8 (e.g., Hopkins & Beacom 2006, Bouwens et al. 
2010). The build-up of stellar mass follows this evolution, as 
does the temporal integral (for z<2, e.g. Ilbert et al. 2010). 
The redshift range z ~ 1 . 5 to 3 has been called the ' epoch of 
galaxy assembly', when about half the stars in the Universe 
form in spheroidal galaxies. It has been shown that star for- 
mation shifts systematically from lower luminosity galaxies 
(Lpjfi ~ 10 10 L ) at low redshift, to high star formation 



rate galaxies (Lfir > 10 11 L Q ) at z ~ 2 ('downsizing', 
e.g., Le Floc'h et al. 2005, Smolcic et al. 2009). Another key 
finding is the tight relation between the star formation rate 
and the stellar mass for star forming galaxies up to redshifts 
z~3 (e.g. Noeske et al. 2007; Elbaz et al. 2007; Daddi et al. 
2007; Magdis et al. 2010). High-redshift starforming galax- 
ies are typically selected through optical/IR studies (e.g., 
Steidel et al. 1999, Daddi et al. 2004, van Dokkum et al. 
2006) 

Besides the star forming galaxies, a population of pas- 
sively evolving, relatively massive (stellar masses ~ 10 10 to 
10 11 Mq) galaxies ('red and dead'), have been detected at 
z > 1, comprising roughly 50% of the galaxies selected in 
near-IR surveys (see, e.g., review by Renzini 2006). These 
galaxies must form the majority of their stars quickly at even 
earlier epochs (e.g., Wiklind et al. 2008). 

Lastly, substantial populations of supermassive black holes| 
and galaxies are now being detected routinely back to first 
light and cosmic reionization (z ~ 6.5). Populations in- 
clude normal star forming galaxies, such as the Ly-a se- 
lected galaxies, with star formation rates of order 10M© 
year -1 (Taniguchi et al. 2005, Ouchi et al. 2009) as well as 
quasar host galaxies, often with star formation rates that are 
100 times higher (e.g., Wang et al 2008). The latter likely 
represent a major star formation epoch for very massive 
galaxies within 1 Gyr of the Big Bang. Most recently, GRBs 
are showing great potential to probe galaxy formation out to 
z > 8 (Tanvir et al. 2009, Salvaterra et al. 2009). 

These observations have led to a new model for galaxy 
formation. As opposed to either cooling of virialized, hot 
halo gas, or major, gas-rich mergers, the dominant mode of 
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star formation during the epoch of galaxy assembly may be 
driven by cold mode accretion, or stream-fed galaxy for- 
mation (e.g., Keres et al. 2005, Dekel et al. 2009). Simu- 
lations suggest that gas accretion in early massive galaxies 
occurs along cold streams from the filamentary intergalac- 
tic medium (IGM) that never shock-heat, but streams onto 
the galaxy at close to the free-fall time. This cool gas forms 
a thick, turbulent rotating disk which very efficiently forms 
stars in a few giant regions in the disk (e.g., Bournaud et al. 
2007). The star forming regions then migrate to the galaxy 
center via dynamical friction and viscosity, forming com- 
pact stellar bulges. The process leads to relatively steady, 
active (~ 100 Mq year -1 ) star formation in galaxies over 
timescales of order 1 Gyr. Subsequent dry mergers at z < 2 
lead to continued total mass build up, and morphological 
evolution, but little subsequent star formation (e.g., Skelton 
et al. 2009, Van der Wei et al. 2009, Robaina et al. 2010). 



2 Molecular gas: the key to testing galaxy 
formation models 

While progress has been impressive, the model above is 
based almost exclusively on optical and near-IR studies of 
the stars, star formation, and ionized gas. There remains a 
major gap in our understanding of galaxy formation, namely,| 
observations of the cool gas, the fuel for star formation in 
galaxies. In essence, current studies probe the products of 
the process of galaxy formation, but miss the source. Nu- 
merous observational and theoretical papers have pointed 
out the crucial need for observations of the cool molecular 
gas feeding star formation in galaxies (e.g., Dressier et al. 
2009; Obreschkow & Rawlings 2009). 

The molecular gas density history of the Universe: 

The key to future studies of galaxy formation is to determine 
the evolution of the molecular gas density of the Universe. 
Over the last two decades, the star formation history of the 
Universe (SFHU) plot has been perhaps the most fundamen- 
tal tool for studying galaxy formation. At the same time, 
very detailed studies of star formation in nearby galaxies 
have reached two critical conclusions. First, star formation 
relates closely with the molecular gas content, but has lit- 
tle relation to the atomic neutral gas content (e.g., Bigiel et 
al. 2009). This fact is verified at high z by the lack of evo- 
lution of the cosmic HI mass density, as determined from 
damped Lya systems, over the same redshift range where 
the star formation rate density is increasing by more than 
an order of magnitude (Prochaska & Wolfe 2009). Second, 
once molecular gas forms, it forms stars, to first order, ac- 
cording to a star formation law (e.g., Kennicutt 1998, Bigiel 
et al. 2009, Leroy et al. 2009, Krumholz et al. 2009, Daddi 
et al. 2010b, Genzel et al. 2010). In low-density environem- 
nts, the general idea is that, once a giant molecular cloud 
becomes self-gravitating, star formation proceeds via local 
processes inherent to the GMC. If the universal star forma- 
tion law were to hold to high redshift, then the classic SFHU 



plot would just be a reflection of the molecular gas density 
history of the Universe. 

The stellar to gas mass density ratio: A related, and 
similarly fundamental issue is to compare the stellar and gas 
mass of galaxies versus redshift. Two extremes scenarios 
can be envisioned. At one extreme, the gas builds up to large 
values (approaching 10 11 M©), and then a dramatic star- 
burst is triggered via, e.g., a major gas rich merger, when the 
gas is rapidly converted into stars on dynamical timescales 
~ 10 8 years. At the other extreme, a continuous supply of 
gas comes from the IGM, cooling into molecular clouds and 
forming stars at a rate comparable to the free-fall rate of gas 
onto the galaxy over ~ Gyr timescales. This latter case cor- 
responds to the cold mode accretion model. This point is 
currently addressed by a number of observational programs 
that aim at studying star forming galaxies at z~ 2 (Daddi et 
al. 2010a, 2010b, Tacconi et al. 2010, Genzel et al. 2010, as 
discussed in Sec. 4.3 and 4.4 below). 

More generally, the study of the cool molecular gas phase| 
in 'typical' high high-redshift galaxies (with SFR^ 100 M yr _1 )| 
provide key insight into the process of galaxy formation in 
a number of ways: 

- A universal star formation law: The low order transi- 
tions of CO provide the most direct means of measur- 
ing the total molecular gas mass. These measurements 
can then be compared to star formation rates to deter- 
mine the evolution of the fundamental star formation 
laws relating cool gas and star formation. Is there a uni- 
versal relationship that governs the efficiency of con- 
version of molecular gas to stars in galaxies of various 
types through cosmic time? What are the fundamental 
physical parameters behind the relationship? For a first 
attempt to address this topic see Sections 4.3 and 4.4 
below. 

- How does the ratio of molecular gas to stellar mass evolve ?| 
This ratio is typically < 0.1 in the nearby Universe, even 
for gas rich galaxies such as the Milky Way. Observa- 
tions of a few the sBzK, and BX/BM galaxies provide 
evidence that this ratio may increase to > 1 at z <~ 2 
(Daddi et al. 2010a, Tacconi et al. 2010). 

- Galaxy dynamics and star formation: Molecular line ob- 
servations are the most direct means of studying galaxy 
dynamics. Are these systems dominated by rotating disks,| 
as predicted by cold mode accretion, or major mergers? 
What is the relative magnitude of ordered motion vs. 
turbulence? What are the stability criteria, and how do 
they compare to the distribution of star formation? 

- Black hole - bulge mass relation: Establishing a galaxy 's| 
mass through molecular gas dynamics can be used to 
test the evolution of the black hole - bulge mass rela- 
tion back to the first galaxies. Indeed, in many cases, 
the molecular line observations are the only means to get 
galaxy dynamics for galaxies with optically bright AGN 
(e.g., Walter et al. 2003, Shields et al. 2006, Riechers 
et al. 2008, 2009). The local relation suggests that the 
formation of black holes and galaxies are closely tied, 
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ular gas in these objects (> 10 10 M©, i.e. comparable to the 
submillimeter galaxies, e.g., Solomon & Vanden Bout 2005, 
Coppin et al. 2008). Quite remarkably, such reservoirs have 
now been detected in a sizable sample of quasars at redshift 
6 and greater, i.e. at the end of cosmic reionization when 



but the origin of this relationship remains a mystery (see 
Sec. 4.2). 

ISM physics and chemistry: Detailed studies of the phys-| 
ical conditions in the ISM of primeval galaxies can be 
performed using numerous molecular line transitions, 
including excitation studies of CO, observations of dense| the Universe was less than one Gyr old (Walter et al. 2003, 



gas tracers like CN, HCN and HCO + , and other key 
astrochemical tracers. Such observations determine the 
temperature and density of the ISM, molecular abun- 
dances and the fraction of dense gas. To date, only a 
few extreme high-redshift systems have been detected 
in these dense gas tracers (e.g. Vanden Bout et al. 2003, 
Riechers et al. 2006, 2007a, 2007b, Garcia-Burillo et al. 
2006) 

3 Recent progress on molecular lines 
observations in early galaxies 

3.1 Submillimeter Galaxies (SMGs) 



2004, Bertoldi et al. 2003, Maiolino et al. 2007, Carilli et al. 
2007, Wang et al. 2010). Resolved imaging of a few quasar 
hosts (Walter et al. 2004, Riechers et al. 2008a, b, Riech- 
ers et al. 2009) imply that the molecular gas is extended on 
kpc scales around the central black hole. It remains puzzling 
how such amounts of cold molecular gas, rich in C and O, 
can be distributed over many kpc scales on short (few 100 
Myr) timescales. 

The high FIR luminosity of the quasars, in conjunc- 
tion with their massive molecular gas reservoirs imply high 
star formation rates in the host galaxies of SFR > 1000 
M© year -1 . This finding is corroborated by the detection of 
[CII] emission in a quasar host at z=6.42 (Maiolino et al. 

2005, Walter et al. 2009). In fact, spatially resolved mea- 
Surveys at submm wavelengths have revealed highly ob- surements of the [CII] emission were able to constrain the 
scured, extreme starburst galaxies (SFR > 1000 M year - *)| size of the starburst to be emerging from the central ~1 .5 kpc| 
at high redshift, the so-called submillimeter galaxies (SMGs,| of the quasar host. This implies extremely high star forma- 
e.g., Smail et al. 2007, Blain et al. 2002). Subsequent stud- tion rate surface densities (~ 1000 M Q yr -1 kpc-2), at the 



ies have detected large reservoirs of molecular gas in these 
systems (> 10 10 M Q ; e.g., Frayer et al. 1998, 1999, Gen- 
zel et al. 2003, Greve et al. 2005, Tacconi et al. 2006, 2008, 
Knudsen et al. 2009, Ivison et al. 2010). General properties 
of the SMGs are summarized in the review by Solomon & 
Vanden Bout (2005) and the detections to date are summa- 
rized in Figure 1 (as discussed in Sec. 4.3). 



maximum what is allowed theoretically based on Edington- 
limited star formation of a radiation-pressure-supported star-| 
burst (Walter et al. 2009a). 

The central black hole masses of the quasars that are 
currently detectable reach masses of many 10 9 M . If the 
local relation seen between the central black hole and sur- 
rounding bulge mass (e.g., Magorrian et al. 1998, Ferrarese 



High-resolution imaging of SMGs by Tacconi et al. (2006,f t al. 2000, Gebhardt et al. 2000) were to hold at these high 



2008) have shown that their molecular gas reservoirs are 
compact, with a median diameter of ^4 kpc. This morphol- 
ogy can be explained if these galaxies are the results of 
mergers, in which the molecular gas settles in the centre 
of two interacting galaxies (leading to a starburst). Recent 
CO imaging of one of the most distant SMGs suggests that 
cold mode accretion may also play a role in powering the 
ongoing star formation (Carilli et al. 2010). 

Chapman et al. (2005) have shown that the median red- 
shift of SMGs is z^2.3 — the existence of a substantial 
high-redshift (z>4) tail of the submillimeter galaxy popu- 
lation that host molecular gas has recently been established 
by a number of studies (Schinnerer et al. 2009, Daddi et al. 
2009a, 2009b, Carilli et al. 2010, Coppin et al. 2010). 

While the study of these SMGs has opened a new win- 
dow on the optically obscured Universe at high redshift, the 
fact remains that these sources are rare, likely represent- 
ing the formation of the most massive galaxies in merger- 
driven hyper-starbursts at high redshift, with gas depletion 
timescales < O.lGyr, 

3.2 Quasars 

Likewise, surveys of molecular gas in quasar host galaxies 
have revealed the presence of massive reservoirs of molec- 



redshifts, one would expect bulge masses of order 10 12 M Q . 
Such high masses exceed dynamical mass estimates of a few 
sources by more than an order of magnitude (Walter et al. 
2004, Riechers et al. 2008a,b, 2009), which in turn indicates 
that black holes built up their masses more quickly than the 
stellar bulges (at least in some quasars). 

3.3 'Normal' starforming galaxies 

More recently, major progress has been made on detecting 
CO emission from more normal star forming galaxies at 
high z. These come in the form of near-IR and color selected 
galaxy samples at z ~ 1.5 to 3. The BzK color selection 
technique in particular has revealed a sample of star forming 
galaxies at z ~ 1.5 to 2.5, selected at near-IR wavelengths, 
with stellar masses in the range 10 10 to 10 11 M , and star 
formation rates ~ 100 M yr -1 (Daddi et al. 2004). These 
galaxies have a space density > 10 times those the submm 
galaxies. 

Daddi et al. (2008) showed that star forming BzK galax- 
ies (sBzK) uniformly contain molecular gas reservoirs of 
comparable mass to the submm galaxies, and yet they are 
forming stars at ~10 times lower rates. This leads to high 
gas depletion times of up to ~1 Gyr (see Fig. 1). Follow- 
up observations showed that massive gas reservoirs have 
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Fig. 1 Comparison of molecular gas masses and total IR bolometric luminosities (taken from Daddi et al. 2010b) for 
all high-redshifted systems detected in CO to date: BzK galaxies (red filled circles; Daddi et al. 2010a), z ~ 0.5 disk 
galaxies (red filled triangles; Salmi et al. in prep.), z = 1-2.3 normal galaxies (Tacconi et al. 2010; brown crosses), SMGs 
(blue empty squares; Greve et al. 2005, Frayer et al. 2005, Daddi et al. 2009a, 2009b), QSOs (green triangles; Riechers et 
al. 2006, Solomon & Vanden Bout 2005), local ULIRGs (black crosses; Solomon et al. 1997), local spirals (black filled 
square: Leroy et al. 2009; black filled triangles: Wilson et al. 2009). The two nearby starbursts M82 and the nucleus of 
NGC 253 are also shown (see Daddi et al. 2010b for references). The solid line (slope of 1.31 in the left panel) is a fit to 
local spirals and BzK galaxies and the dotted line is the same relation shifted in normalization by 1 . 1 dex. For guidance, 
two vertical lines indicate SFR = 2 and 200 M yr _1 in the right panel. For more details see Daddi et al. (2010b). 



been detected in each BzK galaxy that was targeted (Daddi 
et al. 2010a). The implied gas fractions are very high and 
the gas mass in these galaxies is comparable to or larger 
than the stellar mass, and the gas accounts for 50-65% of 
the baryons with the galaxies' half-light radius. New ob- 
servations of an even larger sample of 'normal' starforming 
galaxies have confirmed and extended these findings (Tac- 
coni et al. 2010). 

In a few cases, the molecular gas emission in these ob- 
jects could even be spatially resolved (see Tacconi et al. 
2010, Fig. 2, for a particularly striking example) allowing 
for a determination of the dynamical mass. A good con- 
straint on the dynamical mass also allows one to derive an 
independent measurement of the CO-to-H2 conversion fac- 
tor (as done in the BzK sample by Daddi et al. 2010a). The 
interesting result of such an analysis is that the BzK's have 
a conversion factor similar to the Galaxy, i.e. a factor of ~ 5 
larger than what is found in local ULIRGs. 

In conclusion these systems likely represent the major 
star formation epoch in early massive galaxies (presumably 
driven by cold mode accretion). In the current picture these 
galaxies may turn into passively evolving galaxies at z ~ 1, 
which may eventually turn into 'red and dead' cluster el- 
liptical galaxies seen today. They are substantially differ- 
ent from the quasar and submillimeter population discussed 
above in terms of their larger spatial sizes, lower star for- 



mation efficiencies, and lower molecular gas excitation (e.g. 
Dannerbauer et al. 2009, Aravena et al. 2010). 

3.4 Different Star Formation Laws at High Redshift? 

The detection of molecular gas in systems spanning a wide 
range of parameters now allows one to investigate the star 
formation law (i.e. the dependence of the star formation rate 
on the available gas) at high redshift. In Figure 1 (left), we 
show the summary plot of Daddi et al. 2010b. Here the in- 
frared luminosities (a proxy for the SFR) is plotted as func- 
tion of molecular gas mass (Mh2) of all high-redshift ob- 
jects detected in CO to date, including submillimeter galax- 
ies, quasars and BzK,BM/BX galaxies (see figure caption 
for references). As discussed in Daddi et al. (2010b), two 
different conversion factors (to derive Mh2) have been used 
for the high redhift galaxy populations in this plot: the Galac-| 
tic one for the local spiral galaxies and the BzK/normal 
galaxies at z~ 1.5-2. 5 and a lower 'ULIRG' conversion fac- 
tor for the ULIRS, submillimeter galaxies and quasars. Plot- 
ted this way, there are two sequences emerging (each of 
which can be fitted with a power law of slope ~1.3), one 
for the 'starbursts' (ULRIGS, SMGs and QSOs) and one 
for the 'disk' galaxies. It should be stressed that the off- 
set in sequences seen here is not only due to the choice of 
conversion factors (e.g., Fig. 13 in Daddi et al. 2010a). The 
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right hand panel shows the same data, but the y-axis shows 
the ratio of infrared luminosity to molecular gas mass (i.e. 
the the gas depletion time) plotted as a function of total in- 
frared luminosity. It is obvious from this figure that the gas 
depletion times in the BzKs/normal galaxies are about an 
order of magnitude longer than in the systems that are pre- 
sumably undergoing mergers (such as the ULIRGs, SMGs 
and QSOs; see also discussion in Genzel et al. 2010). These 
different star formation efficiencies appear to be regulated 
by the dynamical timescales of the systems (Daddi et al. 
2010b, Genzel etal. 2010). 

These initial molecular gas observations of the differ- 
ent population of star forming galaxies during the epoch of 
galaxy assembly present a remarkable opportunity to study 
the formation of normal galaxies at critical early epochs, 
when most of their stars form, and when the dominant baryon| 
mass constituent was gas, rather than stars. The galaxies are 
typically extended on scales ~ 1" in the optical, providing 
ideal targets for follow-up high resolution imaging of the 
gas using the EVLA and ALMA. 

4 Outlook: EVLA and ALMA 

Major progress in the studies of the molecular gas content 
of galaxies throughout cosmic times is currently limited by 
the following technological restrictions: 

- Sensitivity: Although remarkable progress has been made| 
in recent years, the sensitivities of current interferome- 
ters do only allow to study the molecular gas content in 
objects that have high star formation rates and/or high 
molecular gas masses (of order 10 10 M Q , SFR^100M Q y 
To push this to more typical object at high redshift (with 
SFR that are an order of magnitude lower) much higher 
sensitivites are required. ALMA, with its unprecedented 
collecting area, will without doubt revolutionize the field| 
of quantifying the molecular gas content in high-z sys- 
tems, e.g. through dedicated studies of molecular line 
deep fields. 

- Resolution: Spatially resolved imaging of galaxies re- 
mains an important task to derive the extent of the molec-| 
ular gas emission and to constrain the dynamical masses 
of the systems (which in turn give independent estimates 
of the conversion factor to derive molecular gas masses 
from CO observations), and to search for observational 
evidence for interactions and/or cold mode accretion. 
Both EVLA and ALMA will enable such studies at high 
spatial resolution - but even with these new facilities 
such studies will be restricted to small samples due to 
the intrinsic faintness of the sources. 

- Spectral Coverage: As galaxies are observed at higher 
and higher redshifts, their molecular lines are redshifted. 
Since the molecular gas excitation has proven to be dif- 
ferent in quasars, submillimeter galaxies and more nor- 
mal galaxies (e.g., Weiss et al. 2005, 2007a, 2007b, Dan- 
nerbauer et al. 2009), it is important that the same trani- 
sions (e.g. rotational J transitions of carbon monoxide) 



are observed to compare the various galaxy populations. 
The EVLA will enable detailed studies of the ground 
transitions of CO up to high redshift. On the other hand 
ALMA will enable studies of molecular gas tracers that 
are currently very difficult to observe (e.g. [CII], [Nil], 
[OIII] and other fine structure lines, e.g., Walter et al. 
2009b). 
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